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ABSTRACT 

I present the results of the search for an optical precursor to the naked-eye burst ~ GRB080319B, which reached 5.87™ optical peak 
luminosity in the "Pi of the Sky" data. A burst of such a high brightness could have been preceded by an optical precursor luminous 
enough to be in detection range of our experiment. The "Pi of the Sky" cameras observed the coordinates of the GRB for about 20 
minutes prior to the explosion, thus provided crucial data for the precursor search. No signal within 3ct limit was found. A limit of 12" 
(V-band equivalent) was set based on the data combined from two cameras, the most robust limit to my knowledge for this precursor. 
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Observations of the optical emission accompanying gamma-ray 
bursts have played a crucial role in understanding this phe- 
nomenon. The measurement of the distance to these explosions 
has allowed us to estimate their true energy and thus significantly 



constrain the models describing their origin (van Paradijs et al. 
IWT 



The presence of very bright optical counterparts during 
the gamma-ray emission was first established by A kerlof et al.| 
1999|l and the existence of prompt optical emission associated 



with the prompt gamma-ray was established in two papers by 



Vestrand et al. (2005, 2006). Optical observations of the naked- 



eye burst GRB080319B shed new light on the prompt emission 
mechanism, revealing for the first time, that a very bright optical 
counterpart became detectable simultaneously (with +5 s preci- 
sion) with the y-ray emission (Racusin et al. 2008||. However, 
much remains unclear about this phenomenon including the na- 
ture of the central engine or interpretation of some of the ob- 
served features, such as y-ray precursors visible for about 20% 
of GRBs ( Lazzati 2005 1. These precursors, as well as prompt 
emission, may coexist with optical precursors. 

The aforementioned precursors have not yet been observed 
and had hardly any chance to be seen due to minimal amount 
of optical observations preceding the GRB explosion. However, 
there are some models predicting the existence of a precur- 
sor. According to these models, precursor detections could ex- 
plain some emission features not yet understood ( jLipunov &| 
|Gorbovskoy 2007) or, more importantly, reveal crucial informa- 



tion about the nature of the GRB central engine (Ouyed et al. 
2002) 1 and answer some long-standing questions. Therefore, the 



discovery of an optical precursor could revolutionise our knowl- 
edge of GRBs. 

Investigations to date of possible optical precursor mecha- 
nisms have reached no firm conclusion. However, in my opin- 
ion, two main types of optical precursors models can be distin- 
guished. The first attributes the precursor generation to the ac- 
tivity of the central engine, such as that produced by a two-stage 
collapse. In one group of these models, the burst is preceded by 
a progenitor star ending its life in a supernova Ibc explosion (jLij 



120081. The explosion should be visible as an optical precursor, 
although its luminosity may be similar to those of standard su- 
peiTiovae and too low, owing to the distance of most of bursts 
observed by optical observatories ( Kuznetsova et al.|2008j ). The 
time between the precursor and the main burst depends strongly 
on the scenario of the explosion, in the case of the SupraNova 
model ( [Vietri &"Stella 199 8) even being as long as months to 
years, but some argue that it can be much shoiter (Dermer 2005J. 
The SupraNova models also allow for optical precursors caused 
by plerionic activity, which can reach 17™ for z ~ I in R-band 



and exist for very short times before the GRB ( Inoue et al. 2003 i. 

Somewhat similar groups of models predict a supernova ex- 
plosion and the creation of a neutron star, which then undergoes 
a "Quark-nova" explosion leading to a GRB (Ouyed et al. 2007 1. 
These models allow for precursor emission seconds before the 
main outburst ( jOuyed et al.||2002| in optical. X-ray, and y-ray 
regimes and the observation of an optical precursor could yield 
crucial information on the delay between the SN and QN. A dif- 
ferent model ( jLipunov & Go rbovsko y[2007) l, still concentrating 
on the activity of the central engine, analyses a collapse of a 
so-called "spinar" - a rotating object for which equilibrium is 
the result of an interaction between centrifugal and gravitational 
forces. The model also permits optical precursors to appear dur- 
ing the lifetimes of observed y-ray precursors. 

The other branch of models predicting an optical precursor 
focuses on the analysis of the outflow itself, not the central en- 
gine. Some models, such as that of Lyutikov & Usov (2000i, 
predict y-ray precursors a with thermal spectrum. However, the 
precursor itself would be 10-100 times less intense than the main 
burst, so its optical tail would be very weak. On the other hand, 
Daigne & Mochkovitch ( 2002[ l found that there may be a non- 
thermal optical precursor, which may however be very weak. 

The detection of GRB080319B enhanced our belief in an 
observable optical precursor, as for the first time this object pro- 
vided the astronomical community with an indisputable proof 
that the initial y emission of the gamma-ray bursts can be ac- 
companied by strong optical emission. The burst was moderate 
in terms of gamma radiation, but extremely luminous, reaching 
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5.3™ - 5.87™ (depending on the experiment) at visible wave- 
lengths. When related to the distance z - 0.937, the luminosity 
is the highest observed to date, beating the previous record set 
by SN 2005ap by a factor of 2.5 ■ 10^. The burst was also excep- 
tionally bright at X-rays. Nevertheless, an extrapolation of y-ray 
spectrum to optical energies leads to fluxes much less intense 
than the observed optical radiation, which rules out the possi- 
bility, that the optical flash is the low energetic tail of the high 
energy emission. This observation together with a possible cor- 
relation between peaks in both bands provides support to calls 
for a revision of prompt emission models ( Racusin et al.|2008| 
[Beskin et alTlorO). 

The predictions of optical precursor models are limited in 
scope, most likely owing to the very few observations provid- 
ing any limits on this type of emission, such as Blake & Bloom 
( |2004| l. Perhaps the main motivation for the search for this pre- 
cursor is the knowledge acquired from studies of the naked-eye 
burst that the optical emission can be simultaneous with the y- 
ray emission. If this is due to the coexistence of emission phe- 
nomena in both bands ( |Kumar & Panaitescu||2008l , there is a 
possibility, that the observed y-ray precursors are also accom- 
panied by optical precursors. In most cases, the observed opti- 
cal emission was less luminous than the y-ray emission and, the 
optical precursor, if accompanying y-ray would be too weak to 
be detected. However, the extreme brightness of GRB080319B 
in the optical band, relative to the moderate intensity in y band 
hints that it could be possible to have an optical precursor bright 
enough to be seen, even if the y-ray precursor was not visible to 
experiments. 

Only two constrains on the optical precursor emission to 
GRB080319B have so far been published. The coordinates of 
the burst were observed for 40 minutes starting more than 2 
hours before the detection by 152 cm telescope located in Loiano 
(Greco et al. 2008). The limiting magnitude of four coadded 10 
minute exposures in filter R is 20.3™. However, this observation 
does not exclude the existence of much brighter, but very short 
precursors. The all-sky monitor RAPTOR-Q has set brightness 
limits in the time period covering known y-ray precursor times 
- 30 minutes prior to the main outburst - to 0.5 Jy (~ 9.63™ in 
V filter). The TORTORA experiment set Umits of 8 5™ fo r 100 
seconds prior to the SWIFT trigger ( Karpov et al.||20(j8 i. The 
presented results improve these limits significantly. 

In all the computations in this article, the Acdm model has 
been used, with Hq = 71, - 0.27, and Qa - 0.73. 
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2. The "Pi of the Sky" project 

While the lack of discovery of an optical precursor to date may 
be attributed to the nature of a GRB itself, the poor constrains 
on the phenomenon occurring between second and dozens of 
minutes before GRB is a consequence of present optical obser- 
vation strategies. Nearly all of the optical GRB observatories are 
follow-up experiments. Detectors rotate to the coordinates of the 
burst after receiving a signal from a GRB orbital experiment, 
such as SWIFT BAT or Fermi GBM. This strategy in addition 
to the small field of view of these experiments make observa- 
tion of the optical prompt emission difficult, being possible only 
for very few long bursts and with very fast devices, and nearly 
rule out the possibility of observing an optical precursor This is 
one of the main reasons why "Pi of the Sky" project adopts a 
completely different strategy. 

The "Pi of the Sky" experiment is designed for continuous 
monitoring of a large fraction of the sky with high time resolu- 
tion (~ 10 s). A real-time analysis of the data stream, based on 



Fig. 1. Comparison of the real (left) and simulated (right) 
GRB080319B images (in peak brightness), the second of which 
was close to the frame corner, where the PSF is most deformed. 
The simulated PSF, used for the precursor search described in 
this paper, reproduces the real PSF in great detail. 



a multi-level triggering system, allows discoveries of GRB op- 
tical counterparts independently of satellite experiments (Burd 
et al.|2005[|Malek et al.|2010[ ). This approach resulted in the au- 
tonomous detection of the naked-eye burst GRB080319B at its 
very beginning ( Racusin et al.|[2008) l. The strategy also allows 
the search for optical precursors to GRBs, similar to those de- 
scribed in this article. 

The GRB080319B was detected by the prototype system 
during its first work period (2004-2009) in Las Campanas 
Observatory in Chile. The prototype consists of two cameras 
placed on a paralactic mount, observing the same part of the sky 
to allow the elimination of optical flashes due to cosmic radia- 
tion hitting the CCD sensor. 

To meet the requirement for monitoring a large fraction of 
the sky, the "Pi of the Sky" apparatus makes use of cameras with 
a very wide field of view - about 20° x 20° each. For stars posi- 
tioned far from the optical axis, this causes significant deforma- 
tions of images, which are much larger than in other astronom- 
ical experiments. This was also the case for GRB080319B, for 
which the position of the burst was in the corner of the frame 
up to fo + 36 s. The observed precursor would therefore also be 
deformed and thus large uncertainties would be introduced into 
standard photometric and signal-searching algorithms. 

To improve the brightness and coordinate measurements, a 
model of the PSF in the "Pi of the Sky" system was created. 
A standard method used to determine the PSF shape is based on 
the superposition of real sky star images. However, in the case of 
very deformed stars, it introduces large uncertainties because the 
exact shape is unknown and the star position is thus uncertain. 

Our approach was therefore different. First, the PSF pro- 
files of a point-like source were precisely measured across the 
CCD chip in laboratory conditions and reconstructed with high 
resolution. Second, a model based on Zernike polynomials was 
fitted to high-resolution profiles (Piotrowski 201 1). Finally, the 
laboratory-based model parameters were refitted to real sky im- 
ages, to accurately reproduce the deformations that are charac- 
teristic of a given camera. 

The final model can be used to improve photometric and as- 
trometric measurements, and fitted to the given coordinates to 
perform a detailed search for a signal. Additionally, the model 
was used to create a sophisticated frame simulator, namely one 
able to reproduce the frame with the expected star PSFs, sev- 
eral types of noises and fluctuations, etc. The simulator can 
be used to verify analysis me thods and design future hardware 
( Piotrowski & Zarnecki|20lT i. 
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Fig. 2. Signal value i at the GRB080319B position relative to 
the nearby reference star, obtained from the PSF profile fit to the 
data from the k2a camera of the "Pi of the Sky" prototype as a 
function of time before the burst t - to (vertical bars). We also 
indicate the limits on the precursor luminosity calculated assum- 
ing no signal (thick lines) and taking the actual signal measured 
for each frame (thin lines). Limits are calculated at both the So- 
level (solid lines) and 95% confidence level (dashed Unes). 



Simulated PSFs obtained in this way are very close to real 
star images even for the most deformed stars, as can be seen in 
fig.[T] As the model was tuned to the limited pixel range around 
the star centre, it does not fully account for the long tails of the 
PSF or a halo appearing around very bright stars. However, these 
effects have a negligible influence on the position and brightness 
measurements, as a vast majority of the signal is measured in the 
accurately modelled part of PSF. 



3. Search for the GRB080319B optical precursor 

The search was performed by fitting the model PSF at GRB co- 
ordinates to all the frames covering 19 minutes prior to the ex- 
plosion, on two cameras (with the internal names k2a and k2d) 
of the "Pi of the Sky" prototype. Figure |2] shows the fitted sig- 
nal value of f for the k2a camera, for all considered frames. To 
suppress the systematic uncertainty, the fitted signal I is referred 
to the scale Ip of a nea rby 7.99™ star (HIP 70912 in Hipparcos 
Catalogue (|ESA||1997|), given in the VT (visual Tycho) filte|^ 
No signal exceeding the 3cr limit has been found for this cam- 
era. The standard approach in the case of "no signal", at least 
in astronomical observations, is to quote a 3cr limit assuming a 
measured value of signal to be zero. The estimation of all the 
limits is based on an assumption that the measurement error for 
a small signal that would be emitted by the optical precursor is 
similar to the fit error for the sky background at this position. 
The uncertainty cr was calculated from a fit of the normal distri- 
bution to the histogram of signal ratio . The result is consistent 
with distributions of signal in three empty control areas in the 
burst proximity that were considered as a cross-check. In this 
approach, the resulting 3(t limiting magnitude for the k2a cam- 
era in the polynomial photometry is 11.67™. At the 95% confi- 
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' The "Pi of the Sky" prototype cameras were not equipped in any 
filter, except for an IR+UV cut filter. This resulted in a relatively wide 
spectral sensitivity of the detector, with average A ^ 586 nm and a 



spread of about 77 nm, closest to the VT and Johnson V filters i Zarnecki 
et al.|2011|l. 



Fig. 3. As in Fig. [ijbut for the k2d camera of the "Pi of the Sky" 
prototype. 



dence level (- 1.96cr), the limiting magnitude for k2a camera is 
12.13™. 

Figure 3 shows the fitted signal value ^ for the k2d camera, 
for all considered frames. The performance of this camera was a 
little bit worse than the k2a, so the background fluctuations are 
higher Thus, the limits are slightly worse, the resulting 3cr lim- 
iting magnitude, assuming no signal, being 11.39™. In this ap- 
proach, any source brighter then 1 1.85™ is excluded at the 95% 
confidence level (CL). 

The method that we used also allowed us to set limits on 
particular frames based on the fitted signal level in the GRB co- 
ordinates (the fit results changing from frame to frame owing 
to the background fluctuation where we seek the signal). The 
limit was numerically calculated according to the formula given 
in Feldman & Cousins| (1998) , which allows us to extract lim- 
its with a well-defined confidence level (CL) for any value of 
the measured signal, including negative signal fluctuations. The 
3cr limits (99.73% CL) calculated on single frames of the k2a 
camera fluctuate in most cases between 11.5™ and 12.25™. The 
95% confidence level limits (which seem more appropriate for 
a single-frame analysis) fluctuate in most cases between 1 1 .6™ 
and 13™. For the k2d camera, the corresponding limits fluctuate 
between 11™ and 12.3™ (3cr limits) and 11.2™ and 13.1™ (95% 
confidence level limits) in most cases. 

The combined scale for a precursor was then computed as 
the weighted average of scales fitted for both cameras. No signal 
above 3cr level was found in the combined signal distribution. 
However, most limits greatly improve, as shown in fig. |4] where 
limits obtained in this analysis are compared to the measured 
GRB080319B Hghtcurve ( [Racusin et al.|2008| l. The standard 3cr 
limit calculated assuming a zero signal is 12™, which is fainter 
by 0.33™ than the single k2a camera limit, and by 0.61™ than 
the single k2d camera limit. The 95% CL limit also increases by 
0.33™ relative to the k2a limit and 0.62™ relative to that of k2d, 
reaching 12.47™. The 3cr limit based on the measured signals for 
single frames is between 11.5™ and 12.6™ for most frames. The 
95% CL limit for most frames is contained between 11.7™ and 
13.3™. 

In general, the combined limits based on signals measured 
in the two cameras are more stable than limits for the single, 
k2a camera. The measurement error was reduced by nearly 30%, 
from 0.01 1 for k2a to 0.008 for the combined signal. This is the 
main reason for the limit improvement with the two cameras. 
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Fig. 4. Limiting magnitude M ( V-band equivalent) for the opti- 
cal precursor emission from GRB080319B as a function of time 
before the burst t - to. Limits to the precursor luminosity were 
calculated assuming no signal (thick lines) and taking the actual 
signal measured for each frame (thin lines). Limits were cal- 
culated at the 3cr level (solid lines) and 95% confidence level 
(dashed lines) and are based on combined data from two cam- 
eras of the "Pi of the Sky" prototype. The points with error bars 
repres ent the GRB080319B "Pi of the Sky" flux ( jRacusin et al.| 



4. Conclusions 

The coordinates of GRB080319B - the naked-eye burst - were 
observed 19 minutes prior to the explosion. The data have been 
analysed with a dedicated signal search algorithm, based on a so- 
phisticated model of the "Pi of the Sky" PSF. Although no signal 
exceeding 3cr limit has been found for either single camera or the 
combined measurements, new, stronger limits on precursor lumi- 
nosity have been established. The combined 3cr limit, assuming 
no signal is a 12™ apparent magnitude, being just 0.35% of the 
peak "Pi of the Sky" flux, which is a significant improvement 
compared to the limit set by RAPTOR-Q in the similar time pe- 
riod, namely 2.5% of the peak flux (Wozniak et al. 2009). The 
rough estimation gives the absolute magnitude precursor limit of 
-31.9™. 

The observation is consistent with the plerionic precursor 
model, which gives ~ 17™ optical precursor in our spectral 
range for GRB at z ~ 1 ( Inoue et al.||2003 1. For the naked-eye 
burst (z = 0.937), this gives a rough estimate of the 16.8™ pre- 
cursor luminosity (with the same cosmological assumptions as 
above), far beyond the "Pi of the Sky" range. The expected su- 
pernova explosions at the GRB080319B distance would also be 
beyond our range, even if the event were superluminous, as in 
the case of the expected quark-novae ( Ouyed et al.|2002: i. 



On the other hand, if we assume that the optical precursor 
luminosity is proportional to the main peak luminosity, and the 
proportion is similar to that observed for y-ray emission, i. e. 10- 



100 times dimmer than the peak luminosity ( [Lyutikov & Usov 
12000) 



we can rule out the possibility of GRB080319B having 
an optical precursor According to our presented calculations, 
the precursor would have to be at least 300 times dimmer than 
the peak luminosity. Future calculations, which were beyond the 
scope of this paper, for specific models of this burst should be 
performed and compared to the aforementioned limit, to provide 
improved constrains on models parameters. 
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